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ABSTRACT 



I We present a synthetic spectral analysis of Far Ultraviolet Spectroscopic Ex- 

^ I plorer (FUSE) and Hubble Space Telescope/Space Telescope Imaging Spectrograph 

' (HST/STIS) spectra of 5 dwarf novae above and below the period gap during qui- 

• ■ escence. We use our synthetic spectral code, including options for the treatment of 

■ the hydrogen quasi- molecular satellite lines (for low temperature stellar atmospheres), 

OO . NLTE approximation (for high temperature stellar atmospheres), and for one system 

(RU Peg) we model the interstellar medium (ISM) molecular and atomic hydrogen lines. 
In all the systems presented here the FUV flux continuum is due to the WD. These 
^ I spectra also exhibit some broad emission lines. In this work we confirm some of the 

' previous FUV analysis results but we also present new results. For 4 systems we com- 

bine the FUSE and STIS spectra to cover a larger wavelength range and to improve 
the spectral fit. This work is part of our broader HST archival research program, in 
which we aim to provide accurate system parameters for cataclysmic variables above 
and below the period gap by combining FUSE and HST FUV spectra. 

Subject headings: Stars: white Dwarfs, Stars: dwarf novae (SS Aur, EY Cyg, VW Hyi, 
RU Peg, EK TrA). 



^Based on observations made with the NASA-CNES-CSA Far Ultraviolet Spectroscopic Explorer. FUSE is oper- 
ated for NASA by the Johns Hopkins University under NASA contract NAS5-32985 

^Visiting at the Space Telescope Science Institute, Baltimore, MD 21218, godon@stsci.edu 
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Introduction 



1.1. Dwarf Novae in Cataclysmic Variables 



Cataclysmic variables (CVs) are short-period, semi-detached compact binaries in which the 
primary star, a white dwarf (WD), accretes matter and angular inomen tum from the secondary, a 
main-sequence star filling its Roche lobe (the mass donor) (jWarner 1 ll995l l. The matter is transferred 
by means of either an accretion disk around the WD, or an accretion column or curtain - when the 
WD has a strong magnetic field. Ongoing accretion at a low rate (quiescence) can be interrupted 
every few weeks to months by intense accretion (outburst) of days to weeks (a dwarf nova accretion 
event), and every few thousand years by a thermonuclear runaway explosion (TNR - the classical 
nova event, due to the ignition of the accreted layer of hydrogen-rich material). CVs are divided in 
sub-classes according to the duration, occurrence and amplitude of their outbursts. The two main 
types of CVs are Dwarf Nova systems (DNs; weakly- or non-magnetic disk systems found mostly 
in their quiescent state which lasts much longer than their outburst), and nova-like (NL) systems 
which form a less homogeneous class. NLs includes disk systems found mostly in their high state, 
intermediate polars (IPs; with a magnetically truncated inner disk), polars (devoid of disk due to 
the strong magnetic field of the WD), and other systems e.g. that never go into outburst or which 
cannot be classified as DNs (jWarner Ill995l ) (this includes the helium systems - AM CVn - with the 
shortest period). Dwarf nova systems includes the U Gem systems, the SU UMa systems, and the 
Z Cam systems. The U Gem systems are the typical DNs, i.e. those systems exhibiting normal DN 
outbursts; the SU UMas exhibit both normal DN outbursts, and superoutbursts, which are both 
longer in duration and higher in luminosity than normal DN outbursts; and the Z Cam systems 
have standstills where they remain in a stat e of intermediate opti cal brightness for a long time 
(here we adopt the classification according to iRitter &: Kolb I (j2003l )). 



Interestingly enough, while the binary orbital period in CV systems ranges from a fraction of an 
hour (AM CVn systems) to about a day (e.g. GK Per), there is a gap in the orbital period between 
2 and 3 hours where almost no systems are found (hereafter the "period gap"). For example U 
Gem and Z Cam DN systems are found above the period gap, while the SU UMa DN systems are 
found below the period gap. It is not known whether the systems are evolving from a longer period 
to a shorter period (across the gap) or whether the systems above the gap are altogether different 
from the systems below the gap. 

The now widely acc e pted interpretation of the quiescence /outburst cycle is the disk instability 
model (DIM, ICannizzo I (119981 )). It is assumed that during the quiescent phase the matter in the 
disk is cold and neutral and the disk is optically thin because of its low density, while during 
outburst the matter in the disk is ionized and becomes optically thick as the mass accretion rate 
within the disk increa ses. The basic principle of the DIM theory depends heavily on the unknown 
viscosity parameter a (jShakura &: Sunvaev Ill973l ) and on the mass transfer rate during the different 
phases. DNs, unlike other CVs, offer a fairly reliable estimate of their dis tances via the a bsolute 
magnitude at maximum versus orbital period relation for DNs found by IWarner I (jl995l ). This 
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relationship is consis tent with theory JCannizzo I Il998l ). The mass accretion rate within the disk 
has been taken from [Patterson I (|1984| ). which is, however, only a first order estimate. In the last 
decade, the disk mass accretion rate of many systems has been deduced more accurately at given 
epochs of outburst or quiescence using spectral fitting techniques. The accretion rate is usually 
a function of time (especially during the outburst itself; and it is also a function of radius r due 
to wind-outflow from the disk) and, consequently, it is difficult to assess its time-averaged value 
accurately. 



Recent advances in theory (jTownslev &: Bildsten 1 12004 ) have shown that the average mass ac- 
cretion rate onto the WD in DNs can be deduced if one knows the mass of the accreting WD and its 
effective surface temperature during quiescence, therefore providing an additional and independent 
way to assess M (or more precisely M(r = R^,) - the mass accretion rate at one stellar radius R^, 
onto the stellar surface, which might be different than the mass accretion rate in the disk M(r) 
at a radius r if there is an outflow from the disk: dM/dr < 0. ). Consequently, in order to put 
more constraints on the theories we need to know the properties (mainly the temperature and mass 
of the WD) of these systems above as well as below the period gap. There is, however, a critical 
shortage in knowledge of the WD properties (effective temperature T^^, gravity Log{g), projected 
rotational velocity Vrotsin{i), chemical abundances, accretion belts?) in DNs above the period gap. 
Thus, detailed comparisons of accreting WDs above and below the gap cannot be made. 



For systems below the gap, with orbital periods near the period minimum (jWarner I Il995l ). 
the distribution of temperatures are centered at ~15,000K with only a narrow range identified at 
present. This distribution appears to manifest the effect of long term compressional heating at a 
time averag ed accretion rate of 2 x IQ-^^MqJt-^ (Townsley and Bildsten 2002; Sion et al. 2003). 
It appears that WD ^ej^'s for systems above the gap are higher than WD temperatures in systems 
below the gap, due to the systems above the gap having larger disks (with higher mass transfer 
rates) and more massive (somewhat earlier-type) secondaries. Some disks may remain optically 
thick even during quiescence so that the WDs are heated to a greater extent than systems below 
the gap. It is not yet known whether the WDs in systems above the gap are rotating more slowly 
than WDs in systems below the gap where presumably the CVs are older and with a longer history 
of angular momentum transfer via disk accretion. Thus far the only DNs above the gap whose 
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total of 2 Z Cam systems, 4 U Gem systems and the peculiar DN WW Ceti (although classified as 
a U Gem, others suggest it is a Z Cam). 

Under our broader iJST archival research we have started to secure accurate system parameters 
(M, i, M^di '^eff, Vrotsin{i), chemical abundances,etc..) for CVs (DNs and NLs) above and 
below the period gap by fitting synthetic spectral models (WDs and accretion disks) to combined 
FUSE+HST spectra from the MAST archive. We have identified 25 CV systems for which the 
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FUSE and HST (STIS, FOS or GHRS) spectra match and can be combined. As a part of this HST 
archival research, we present in this paper the analysis of HST/STIS and {FUSE) archival spectra 
of 5 DNs during quiescence using synthetic stellar spectra together with the FUSE spectrum of RU 
Peg. 



1.2. Five Dwarf Nova Systems 

The five DNs are listed in Table 1 with their system parameters as follows: column (1) Name, 
(2) CV subtype, (3) reddening value E(B-V), (4) distance, (5) orbital period in days, (6) orbital 
inclination in degrees, (7) spectral type of the secondary, (8) mass of the primary in solar masses, 
(9) mass of the secondary in solar masses, (10) apparent magnitude in outburst, and (11) apparent 
magnitude in quiescence. The references are listed below the table. Three systems are above the 
period gap: FY Cyg, SS Aur and RU Peg; and two systems are below the period gap: VW Hyi 
and FK TrA. 

Among these 5 objects, 4 DNs in quiescence are directly from our HST archival program with 
matching FUSE and HST/STIS spectra: FY CYG, SS Aur, VW Hyi, and FK TrA. We supple- 
ment this study here with the analysis of RU Peg's FUSE spectrum. For SS Aur, VW Hyi and 
FK TrA we assess the WD parameters by modeling for the first time for these 3 objects the com- 
bined FUSE+STIS spectra. For RU Peg's FUSE spectrum, and FY Cyg's combined FUSE+STIS 
spectrum we improve the fittings that were carried out previously. 

For all the spectra we present here, we use the latest versions of the stellar model atmospheres 
and synthetic spectra codes (see section 3); this includes options for the treatment of the hydrogen 
molecular satellite lines (for modeling cooler WDs), and NLTF atmos phere models (fo r mode ling 



hotter WDs). We also adopt the extinction values F(B-V) given in iBruch Sz Fngel I (| 19941 ) for 
the 5 objects and we deredden the spectra accordingly. For some systems we identify the ISM 
molecular absorption lines in the FUSE spectra; this greatly helps improve the fitting of the chemical 
abundances and WD's temperature. RU Peg has stronger ISM absorption lines and, for this system 
only, we model the ISM atomic and molecular hydrogen opacities. For these 5 systems we obtain 
the temperature, projected rotational velocity and chemical abundances (of C, N, S and Si) of the 
exposed WD atmosphere. We also give upper and lower limits for the mass of the WD and the 
distance to the system. 



2. The Archival Spectra 

The observations log is presented in Table 2. From the AAVSO (American Association of 
Variable Stars Observers) data, we found that all the systems were in optical quiescence at the 
time of the observation (see the last column of Table 2). The STIS spectra of FY Cyg and SS Aur 
are snapshots (lasting 600-700s) with a lower resolution and S/N than the STIS spectra of VW Hyi 
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and EK TVA. With fluxes < 1 x 10"^^ergs s"^cm"2A"\ both EY Cyg and EK TVA are actually 
weak FUSE sources, and RU Peg (though with a FUSE flux ten times higher) is under-exposed 
with an exposure time of only 2800s. 



2.1. The FUSE Archival Spectra 

2.1.1. Processing the FUSE Data 

FUSE'S optical system consists of four optical telescopes (mirrors), each connected to a Row- 
land spectrograph. The four diffraction gratings of the four spectrographs produce four independent 
spectra on two detectors. Two mirrors and two gratings are coated with SiC to provide wavelength 
coverage below 1020 A, and the other two mirrors and gratings are coated with Al and LiF. The 
Al-|-LiF coating provides about twice the reflectivity of SiC at wavelengths >1050 A, and very little 
reflectivity below 1020 A. These are known as the SiCl, SiC2, LiFl and LiF2 channels. 

All the FUSE spectra presented here were obtained through the 30"x30" LWRS Large Square 



Aper ture in TIME TAG mode. The data were processed with CalFUSE version 3.0.7 (jPixon et al. 



20071 ). which automatically handles event bursts. Event bursts are short periods during an exposure 
when high count rates are registered on one or more detectors. The bursts exhibit a complex 
pattern on the detector, their cause, however, is yet unknown (it has been confirmed that they are 
not detector effects). The main change from previous versions of CalFUSE is that now the data 
are maintained as a photon list (the intermediate data file - IDF) throughout the pipeline. Bad 
photons are flagged but not discarded, so the user can examine, filter, and combine data without 
re-running the pipeline. A number of design changes enable the new pipeline to run faster and use 
less disk space than before. Processing time with CalFUSE has decreased by a factor of up to 10. 

To process FUSE data, we follow the same pr ocedure used previously for the analysis of other 



systems (such as WW Ceti, iGodon et al. I (j2006bl )): consequently we give only a short account of 



this procedure. The spectral regions covered by the spectral channels overlap, and these overlap 
regions are then used to renormalize the spectra in the SiCl, LiF2, and SiC2 channels to the flux in 
the LiFl channel. We then produced a final spectrum that covers the full FUSE wavelength range 
905 — 1187 A. The low sensitivity portions of each channel were discarded. In most channels there 
exists a narrow dark stripe of decreased flux in the spectra running in the dispersion direction. 
This stripe has been known as the "worm" and it can attenuate as much as 50% of the incident 
light in the affected portions of the spectrum; - this is due to shadows thrown by the wires on 
the grid above the detector. Because of the temporal changes in the strength and position of 
the "worm", CALFUSE cannot correct target fluxes for its presence. Therefore, we carried out a 
visual inspection of the FUSE channels to locate the worm and we manually discarded the portion 
of the spectrum affected by the worm. We combined the individual exposures and channels to 
create a time-averaged spectrum weighting the flux in each output datum by the exposure time 
and sensitivity of the input exposure and channel of origin. 
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2.1.2. The FUSE Lines 

The FUSE spectra of DNs in quiescence exhibit mainly absorption hnes from the WD itself, 
as the exposed WD is the main FUV component of the system. A second FUV component is 
sometimes present as a flat and featureless continuum contributing in the very short wavelengths of 
FUSE (A < 970A ). This second component could be a hot region of the accretion disk (the inner 
disk) or the boundary /spread layer. Sharp absorption lines from circumstellar (or circumbinary) 
material and/or the ISM are also often seen and, if the system is a weak FUSE source, sharp 
emission lines from air glow (geo- and helio-coronal in origin) are present 

Since the WD is the main component of the FUV spectra of DNs in quiescence, the main 
characteristic of the FUSE spectra is the broad Ly/? absorption feature due to the gravity {Log{g) ~ 
8) and temperature (15, 000K< T < 25, OOOK) of the exposed WD (in the disk this feature is usually 
smoothed out due to velocity broadening, unless the disk is almost face-on). The other most 
common absorption features observed in the spectra of DNs WDs are due to Cm (1175 A), Cii 
(1066 A), Sim (f«1108-1114 A and Rill40-1144 A), and Nil (1085 A when not contaminated by air 
glow). At higher temperatures (T> 25, OOOK), as the continuum rises in the shorter wavelengths, 
the higher orders of the Lyman series also become visible; however, they become narrower. At 
these temperatures the SiV (1073 A) absorption line starts to appear, and, as there is more flux 
in the shorter wavelengths, the C II (1010 A) absorption line also becomes visible. At still higher 
temperature (T> 50, OOOK), the C ii and Si iii lines disappear, and the spectrum becomes dominated 
by high order ionization lines such as Niv (ps923 A), Svi (933.5 & 944.5 A), Siv (1063 & 1073 A), 
Siiv (1066.6 A), and Oiv (1067.8 A). 

On top of the spectrum of the WD, broad emission lines are found in quiescent DN systems, 
usually the Ovi doublet and Cm (977A and 1175 A). These, together with broad emission lines 
from Niv {^923 A) and S vi (933.5 & 944.5 A) are often observed in nova-like systems (e.g. such 
as AE Aqr, V347 Pup, DW UMa, see the MAST FUSE archives) and are not usually present in 
low-inclination DN systems in quiescence. This implies that the gas could possibly be heated by a 
shock (e.g. in the boundary layer) and the broadening and variatio n of the lines su g gests they are 



originating in the disk. This scenario has recently been supported bv lKromer et al. I (|2007l ) who for 
the first time modeled the broad H and He emission lines ab initio as irradiation of the inner disk 
by the hot boundary layer and/or white dwarf. One fully expects that the metallic emission lines 
form in the same way. The DN system here that exhibits the strongest emission lines is EK TrA 
(i=58°). RU Peg (i=33°) does have some strong emission lines too, but lacks a STIS spectrum for 
complete comparison; however, its lUE spectrum clearly show strong emission lines such as C IV 
(1550A) and Siiv (1400A). 



The FUSE spectra of the DN systems (and CVs in general) often show some ISM molecular 
hydrogen absorption, which appears as sharp lines at almost equal intervals (12A) starting at 
wavelengths around lllOA and continuing towards shorter wavelengths all the way down to the 
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hydrogen cut-off around 915A . In the affected FUSE spectra, we identified the most prominent 

molecular hydrogen absorption lines by their band (Werner or Lyman), upper vibrational level 
(1-16), and rotational transition (R, P, or Q) with lower rotational state (J=l,2,3). 

2.2. The HST/STIS Archival Spectra 

2.2.1. Processing the HST/STIS Archival Data 

All the STIS spectra were processed with CALSTIS version 2.19. Except for the spectrum 

of EK TrA (obtained in TIME TAG operation mode), all the spectra were obtained in ACCUM 
operation mode. All the STIS spectra used the FUV MAMA detector and all were centered on the 
wavelength 1425A . The STIS snapshots (with an exposure time of about 600-700s, for EY Cyg 
and SS Aur) were obtained through the 52x0.2 aperture using the G140L optical element. These 
snapshots consist of one spectrum. The STIS spectra of VW Hyi and EK TrA (with an exposure 
time of 4,000s and above) were obtained through the 0.2x0.2 aperture using the E140M optical 
elements, and with 42 echellc spectra each, are of a much higher resolution than the snapshots. 
Towards the longer wavelengths, the echelle spectra do not overlap and five gaps are apparent 
around A «1634A, 1653A, 1672A, 169lA, and 1710A. 

2.2.2. The HST/STIS Lines 

The STIS spectra of quiescent DNs are also dominated by the WD and their main characteristic 
is the broad Lyo absorption feature at ~1216 A. The other very common absorption lines are the 
carbon lines Cl (1266 A, 1561 A, 1657 A), Cll (1335 A), and Cm (1175 A); and the silicon lines 
Sill (1260 A, 1300 A, 1530 A), Sim (1300 A), and Siiv (1400 A). 

The HST/STIS spectrum also exhibits some broad emission lines. The most prominent ones 
are Cm (1175A), Civ (1550A), Sim (1206A), Nv (1240A), Hell (1640A Balmer a), and SilV 
(1400A). With their broad emission lines, both the FUSE and HST/STIS spectra show evidence of 
hot gas. 

2.3. Preparation of the Spectra 

For four objects (EY Cyg, SS Aur, VW Hyi and EK TrA) we wish to combine the FUSE 
spectrum with the HST/STIS spectrum in order to have a larger wavelength coverage and improve 
the synthetic spectral modeling. 

However, for each object, before combining the FUSE spectrum with the STIS spectrum we 
must make sure that (i) the system was observed in the same (high or low) state; (ii) the shape of 
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the spectra (continuum and lines) is similar (if the spectra are not too noisy); (iii) the flux level in 
the FUSE and HST/STIS spectra do not differ more than about 50% (FUSE and HST/STIS arc 
very different instruments, we do not expect their flux level to match exactly); (iv) for long period 
systems for which only short exposure time spectra exist, the FUSE and STIS spectra have most 
probably been obtained at a different binary orbital phase and it is not clear one can combine such 
spectra. We address these points here in detail as follows. 

(i) The observations listed here were not coordinated and consequently some were carried out 
during early quiescence while others were carried out in late quiescence. For EY Cyg, the outburst 
recurrence time is believed to be about 2000 days and the FUSE and STIS spectra were obtained 
in deep quiescence. EK TrA has also a long outburst recurrence time and it was observed in deep 
quiescence (45 days and 155 days after outburst). For SS Aur, both spectra were obtained about 
a month into quiescence. For VW Hyi the FUSE spectrum was obtained 11 days into quiescence 
while the main STIS spectrum was obtained 15 days into quiescence (we discuss VW Hyi in some 
more detail in the results section). 

(ii) For all the systems, a visual analysis shows that the shape of the continuum and lines is 
similar in the FUSE and STIS spectra. 

(iii) Only 2 systems (EY Cyg and SS Aur) have FUSE and STIS flux levels that do not match 
exactly. For EY Cyg the STIS spectrum had to be multiplied by 1.5 to match the flux level of the 
continuum in the FUSE spectrum, while for SS Aur it was multiplied by 0.93. At temperatures of 
about 25,000K a change of 2,000K can produce a change of 50 percent in the flux; the accuracy of 
the temperature we obtained for EY Cyg is therefore of this order (i.e. ±2,000K). 

(iv) Since the systems were observed in quiescence, we expect mainly to see the exposed WD 
with basically no contribution or occultation from the accretion disk (none of these systems is 
eclipsing). The spectra, whether observed at a given phase or averaged over an entire orbit, should 
not differ significantly. Only for spectra with a high S/N ratio and a high inclination do we expect 
the absorption and emission lines to exhibit some red- or blue- shift if they were obtained at a given 
orbital phase, while the orbit-averaged spectra will have broader absorption and emission features. 

Except for RU Peg, all the FUSE spectra exposure times cover a significant fraction of the 
orbital period (and in some cases several periods); therefore, the FUSE spectra are averaged over a 
significantly large fraction of the orbit. The STIS spectra of VW Hyi and EK TrA arc also averaged 
over the orbit, and for these two systems one can therefore combine the FUSE and STIS spectra. 

The STIS snapshots of EY Cyg and SS Aur are very short (less than Iksec) and were taken 
at a particular orbital phase. One might therefore argue that the STIS and FUSE spectra of EY 
Cyg and SS Aur cannot be combined because of that. However, for both EY Cyg and SS Aur, they 
fit except that the absorption lines in the STIS range are not as deep as in the FUSE range. The 
effect of observing the systems at a particular orbital phase is not very pronounced, and this might 
be due to the lower resolution of the STIS snapshots (in comparison to VW Hyi and EK TrA's 
STIS spectra) and the moderate inclination of EY Cyg and SS Aur. Therefore, for EY Cyg and 
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SS Aur, we also combine the FUSE and STIS spectra without concern about the orbital phase. 



The only object for which we find a significant shift in the lines is the long period system RU 
Peg, for which we only have a FUSE spectrum of 2,800s. 



Final l y, we deredden the spectra a ccording to the E(B-V) values given in IVerbunt I (|l987l ) 



La Dous I (|l99lh : lBruch Engel I (|l994l ^. 



3. Spectral Modeling 

3.1. The Synthetic Stellar Spectral Codes 

We create model spectra for high-gravity stellar atmospheres using codes TLUSTY and SYN- 
fflubeny 1ll988l :l iHubeny &: Lanz Atmospheric structure was computed (using TLUSTY) 

assuming a H-He LTE atmosphere; the other species were added in the spectrum synthesis stage 
using SYNSPEC. For hot models (say T > 50, OOOK) we switched the approximate NLTE treatment 
option in SYNSPEC (this allows to consider and approximate NLTE treatment even for LTE mod- 
els generated by TLUSTY). We generate photospheric models with effective temperatures ranging 
from 12,000K to 75,000K in increments of about 10 percent (e.g. 1,000K for TRil5,000K and 
5, OOOK for TRi70,000K). We chose values of Log{g) ranging between 7.5 and 9.5 for consistency 
with the observed mass. We also varied the stellar rotational velocity Vrotsin{i) from 100km s^^ 
to 1000km s~^ in steps of 100km s^^ (or smaller if needed). In order to try and fit the absorption 
features of the spectrum, we also vary the chemical abundances of C, N, S and Si. For any WD mass 
there is a corr esponding radius, or equival ently o ne single yalue o f Loq{q) (e.g. see th e mass radius 



(j200nl ) for different 



relation from iHamada k. Salpeter I (Il96lh or see IWood I (jl99d ) ; iPanei et al. 
composition and non-zero temperature WDs). 

Our suite of stellar spectra generator codes has been implemented and includes also the treat- 
ment of the quasi-molecular satellite Lyman lines. The satellite hydrogen lines appear as stron g 
absorption features near 1400 and I6OOA (Lya, in the lUE and STIS range; iKoester et al. I (|l985l ): 
Nelan &: Wegner J 1985 V) and are so mewhat weaker near 10 60 an d 1078A (Ly/?, in the FUSE range; 



Dupuis et a. 



e.g 



by lAllard et al. 



(I2OO6II). Following iDupuis et al. I (120031 . 120061 ) . we used opacity tables computed 



(Il998l . 11994 l2004al ) to take into account the quasi-molecular satellite line opac- 



ities of H2 and for Lya, Ly/3, and Ly7. The quasi-molecular lines are expected to appear 
when the temperature of the WD is ~20,000K or lower, though the gravity, pressure and magnetic 
field of the WD can also play a role in the f o rmati o n of these lines ( l arger electronic density will 
favor recombination, see e.g. 
discussion). 



Dupuis et al. I (j2003l ): iGansicke et al. I (|2006l ) for a more complete 



1 |http://nova.astro.umd!edul TLUSTY version 200, SYNSPEC version 48 
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3.2. Modeling the ISM Hydrogen Absorption Lines 



For all the systems showing ISM atomic and molecular hydrogen absorption lines, we identify 
these lines in the figures to avoid confusing them with the WD lines. For RU Peg, however, some 
ISM lines are deep and broad and we decided to model them, especially since some of the WD lines 
(such as S IV AA1062.6 & 1073) are located at almost the same wavelengths. 

For RU Peg only, we model the ISM hydrogen absorption lines to assess the atomic and 
molecular column densities. This enables us to differentiate between the WD lines and the ISM 
lines, and helps improve the WD spectral fit. The ISM spectra models are generated using a program 
developed by P.E. Barrett. This program uses a custom spectral fitting package to estimate the 
temperature and density of the interstellar absorption lines of atomic and molecular hydrogen. 
The ISM model assumes that the temperature, bulk velocity, and turbulent velocity of the medium 
are the same for all atomic and molecular species, whereas the densities of atomic and molecular 
hydrogen, and the ratios of deuterium to hydrogen and meta ls (includ i ng he l ium) t o hydrogen can 
be adju sted independe n tly. T he model uses atomic data of iMorton I (l2000l . 120031 ) and molecular 



checked against those of iMcCandliss I (120031 ) 



data of lAbgrall et al. I (1200011. The o ptical depth calculations of molecular hydrogen have been 



The ratios of metals to hydrogen and deuterium to hydrogen are fixed at and 2 x 10~^, 
respectively, because of the low signal-to-noise ratio data. The wings of the atomic lines are 
used to estimate the density of atomic hydrogen and the depth of the unsaturated molecular lines 
for molecular hydrogen. The temperature and turbulent velocity of the medium are primarily 
determined from the lines of molecular hydrogen when the ISM temperatures are < 250K. 

The ISM absorption features are best modeled and displayed when the theoretical ISM model 
(transmission values) is combined with a synthetic spectrum for the object (namely a WD synthetic 
spectrum). 



3.3. Synthetic Spectral Model Fitting 

Before carrying out a synthetic spectral fit of the spectra, we masked portions of the spectra 
with strong emission lines, strong ISM molecular absorption lines, detector noise and air glow. 
These regions of the spectra are somewhat different for each object and are not included in the 
fitting. The regions excluded from the fit are in blue in Figures 1, 3, 5, 7, and 10. The excluded 
ISM quasi-molecular absorption lines are marked with vertical labels in figures 2, 4, and 9. For RU 
Peg, we model the ISM absorption features and the WD separately; namely, we also mask the ISM 
absorption features when we model the WD. 



After having generated grids of models for each target, we use FIT ( Press et al.lll992l ). a 
minimization routine, to compute the reduced Xu (x^ P^r number of degrees of freedom) and scale 
factor values for each model fit. While we use a minimization technique, we do not blindly 
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select the least models, but we examine the models that best fit some of the features such as 
absorption lines (see the fit to the FC/S'E' spectrum alone) and, when possible, the slope of the wings 
of the broad Lyman absorption features. We also select the models that are in agreement with the 
known distance of the system. 

However, in the model fitting, for a given WD mass and radius, the resulting temperature 
depends mainly on the shape of the spectrum in the FUSE range. The flux level at lOOOA (between 
Ly(5 and Ly7) is close to zero for temperatures below 18,000K, at 30,000K it is about 50% of the 
continuum level at llOOA and it reaches 100% for T>45,000K. At higher temperature (T>50,000K) 
the spectrum becomes pretty flat and there is not much difference in the shape of the spectrum 
between (say) a 50,000K and a 80,000K model. When fitting the shape of the spectrum in such a 
manner, an accuracy of about 500-1, OOOK is obtained, due to the S/N. In theory, a fine tuning of 
the temperature (say to an accuracy of about zbSOK) can be carried out by fitting the flux levels 
such that the distance to the system (if known) is obtained accurately. However, the fitting to the 
distance depends strongly on the radius (and therefore the mass) of the WD. In all the systems 
presented here the error on the measured mass of the WD is so large (see Table 1) that it produces 
an error in the temperature much larger than 1,000K. This is because the Ly g and Ly/3 profiles 



depend on both temperature and gravity. For example iGansicke et al. I (j200ll ) derived a best fit 
temperature for EK Tra T(,ff ~ 2360 x Log{g) — 95, where g is the (unknown) surface gravity 
of the WD (7.0 < Log{g) < 9.0). Additional errors are further introduced as the distance and 
reddening are rarely known accurately, therefore increasing significantly the inaccuracy of assessing 
the temperature by scaling the synthetic flux to the observed flux. For these, we decided to fit the 
temperature of each model, based on the shape of the FUSE spectrum in the shorter wavelengths, 
for a temperature accuracy of only about 5%, namely accurate to within 500K for T< 20, OOOK; 
1,000K for 20,000K<T<35,000K; and 2,000K for T>40,000K (see the Results section). While the 
error in temperature is only 500K for a 18, OOOK model and 2, OOOK for a 50, OOOK model, the relative 
error AT/T is the same for all the models. 

For all models we first tried fitting solar abundance models and then tried changing the abun- 
dances of some species to improve the spectral fit. In particular, for each system we fit obvious 
absorption features which can determine the abundance of some specific element (C, Si, S, N). For 
example, the Cii (s^^lOGS A ), the Sim (r^IIIO A ) and the Si-C blend (^1138-1146 A ) are the 
main absorption features in the FUSE range of a T=20,000K WD that help assess the Si and the 
C abundances (these are discussed in detail for each system in the next section). 

The WD rotation {Vrotsin{i)) rate is determined by fitting the WD model to the spectrum 
while paying careful attention to the line profiles in the FUSE portion of the combined spectrum. 
We did not carry out separate fits to individual lines but rather tried to fit the lines and continuum 
in the same fit while paying careful attention to the absorption lines. 



It is important to note that, in the modeling, the depth of the absorption features depends 
not only on the abundances but also on the rotational velocity. Increasing the rotational velocity 
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decreases the depth of the absorption features, but broadens the wings - the net effect being that 
higher metal abundances result with faster rotational velocity. With sufficiently high S/N data, it 
becomes possible to assess the rotational velocity and abundances independently. 



4. Results and Discussion 



All the results are presented in Table 5. For some systems for which the mass of the WD is 
relatively unknown, we list more than one optimized fit result, assuming different values of Log{g) 
(i.e. for each value of Log{g) we find an optimized fit). In the first column we list the name of the 
system; in the second column the assumed surface gravity on a logarithmic scale Log{g) (where 
it is understood that Log is the logarithm to the base 10 - as opposed to log which denotes the 
logarithm to the base e) in cgs units; in column (3) we list the effective surface temperature we 
obtained for the WD; in column (4) we list the projected rotational velocity which was obtained 
when matching the stellar absorption lines; in columns (5), (6), (7) and (8) we list the abundances 
of Carbon, Silicon, Nitrogen and Sulfur (respectively) in solar units; in column (9) we list the 
distance we obtained from the fitting; in column (10) we list the parameter by which the STIS flux 
had to be multiplied to match the FUSE flux; in column (11) we give the xt value; and in column 
(12) we indicate the number of the figure displaying the model fit. 



4.1. EY Cygni 



EY Cygni is a peculiar long period U Geminorum-type of DN, with a very massive accreting 
WD (Ml = 1.26Mq or possibly larger), an outburst ampl itude of 4. 1 rnagnitude with a recurrence 



time of 2000 days, and a n orbital period of 1 1.0238 hrs ([Smith et al. 1119971 : ICostero et al. Ill998 



Tovmassian et al. I l2002l : ICostero et al. 1 12004| ) . Ha and S il ground-base d imaging has revealed it 



is associated with a non-homogeneous shell with a size of about 25arcsec ([Tovmassian et al 



2002 



Sion et al. M2004bl ). It is possible that this shell is the result of a recent nova explosion. The 
spectral type of the secondary is not certain, it is believed to be a K5V to MOV star, from which 
the lower limit of the distance to EY Cyg has been inferr ed to be at least 250pc. Si nce EY Cyg 
is not embedded in the background Cygnus superbubble (jBochkarev fc Sitnik I Il985l ). the upper 
limit for its distance is 700pc. An additional ch a racter istic which makes EY Cyg a peculiar DN 
is its anomalous Nv/Civ ratio. IWinter &: Sion I (j200ll ) first reported a very strong Nv emission 
and a very weak (or absent) C iv, which is atypical of most DNs in which C iv is usually the mos t 



prominent and common emission line in their FUV spectra. More recently iGansicke et al. I (|2003l ) 



reported an analysis of HST/STIS snapshots of 4 CVs (including EY Cyg) with anomalously large 
N v/C IV line flux ratios, similar to those observed in AE Aqr. So far 10 systems have been identified 
with such an anomalous N v/C IV ratio. 



EY Cyg is an extremely weak FUSE source and the spectrum is therefore very noisy, which 
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is the reason we binned its FUSE spectrum at O.sA for the fitting (Figure 1). The STIS snapshot 
spectrum was binned at the default value of O.SSA. Consequently the relative weight of the STIS 
spectrum is about twice as large as that of the F?7S'£' spectrum. In preparation for the model fitting 
we had to multiply the STIS flux by 1.5 to match the FUSE flux level; this change of 50% in flux 
corresponds to a (0.5)4 Ri 10% in the value of T. Therefore, the relative error in temperature in 
our spectral modeling of the combined FUSE+STIS spectrum of EY Cyg is of the order of ~ 0.1 
ab initio. Since the distance to the system is relatively unknown, there is no restriction on the flux 
level to minimize the number of best fits to one in the Log[g) — T^ff parameter space. Therefore, 
we restricted the modeling assuming only two different values for the mass of the WD. Because 
of the extremely large WD mass of the system (which within the margin of error well exceeds the 
Chandrasekhar mass limit for a WD) we ran models assuming Log{g) = 9.0 (corresponding to 
My^d = I.2IM0) and Log{g) = 9.5 (corresponding to M^^a = I.35M0). 

Within the range of the temperature error (10%) the best fit models spread between about 
T=27,000K and T=33,000K. However, the high temperature (T > 32, OOOK) models do not produce 
enough flux in the longer wavelengths of STIS, while the lower temperature (T < 28, OOOK) models 
do not have enough flux in the shorter wavelengths of the FUSE portion of the spectrum. Because 
of that we rejected the T = 27, OOOK and T = 33, OOOK models. In addition, in the observed 
STIS spectrum the bottom of the Lya does not go to zero and this corresponds to models with 
T ~ 30, OOOK and higher. Because of the difference in fluxes between STIS and FUSE we decided 
to check how the results would change if we flt the spectra separately. We found that the STIS 
spectrum gave the same temperature, however with a distance about 20 percent larger as the main 
driving elements in choosing the best fit model were the same, namely: the Lya profile and the 
longer wavelengths of STIS. Fitting the FUSE spectrum alone gave a temperature about 2, OOOK 
higher than for the combined spectrum. We therefore inferred for EY Cyg T = 30, 000Kib2, OOOK 
assuming Log{g) = 9.0. The 2, OOOK error in T produces an error of ib60pc in the distance as listed 
in Table 5. This 30,000K model is shown in Figure 1. In Figure 2 the 30,000K WD model assuming 
Log{g) = 9.0 is shown with the FUSE spectrum binned here at 0.1 A to show the sharp emission 
and absorption lines. 

We also checked how the inclusion of the quasi-molecular hydrogen aflFects the result and found 
that it does not improve the fits. 

Next, we assume Log{g) = 9.5 to check the effect of larger WD mass. Using the same consid- 
erations as before we found that the best fit models range between 30, OOOK and 34, OOOK, or about 
2,000K higher than the Log{g) = 9.0 models. 

The projected rotational velocity we obtained for all the models is relatively low: 100km s~^. 
We confirm the anomalously low Carbon abundance, namely we found that the best Carbon line fits 
are obtained for abundances ranging between 0.01 and 0.05 solar. The low carbon abundance was 
set to match the FUSE absorption features of Cm 1175 A, Cii 1066 A, and Ci-Ciii at 1128 and 
1140 A. There is also a complete absence of C lll emission at 977A and II75A, and C IV emission at 
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1550A . There is a strong and broad N v 1240 A emission line. The Nitrogen abundance (2x solar) 
was obtained by fitting the Ni (1135A), Nil (1185A), and Niii (1189A) lines. However, these lines 
are not very pronounced and are afi'ected by air glow emission. The Sulfur abundance (6x solar) 
was obtained by fitting the Siv (1073A) and Si (1097A) lines. The Silicon abundance (0.6x solar) 
was obtained by fitting the Siiv (1125A), and Sim (1140-1145A) lines. 



Sion et al. I (|2004b|) modeled the combined HST/STIS and FUSE spectra of EY Cygni. Despite 
the noisy, weak FUSE spectrum, they managed to be obtain two possible solutions which in a 
statistical {Chi"^) sense differ little from each other. If the combined FUV spectra was due to a single 
temperature white dwarf alone, then for their adopted mass of I.2M0 (assuming Log{g) = 9) with 
no interstellar reddening (see below), then the WD has Tg^jf = 24, OOOK. However, the combined 
FUSE plus STIS spectra also yielded a similar statistical fit if a two-temperature white dwarf 
model was applied to the data. This fit corresponded to a cooler, slowly rotating WD photosphere 
{Teff = 22, OOOK) with the second temperature iTf.ff = 36, OOOK) representing a putative hot, 
rapidly rotating, equatorial accretion zone (an accretion belt). We note here that their noisy FUSE 
spectrum showed a sharp increase in flux extending even beyond the Lyman limit. This extra fiux 
is an artifact of the CaLFUSE background subtraction and, therefore, in the present woork we 
process the FUSE data with latest (and final) version of CalFUSE (v3.2) which lowered the fiux 
at wavelengths j 950A, removing the shortward fiux excess. The new combined HST/STIS plus 
reprocessed FUSE data is fit in the present paper with the FUSE spectral region shortward of 
950Amasked in the fitting. 

Our single WD component best fit, assuming Log{g) = 9.0, has T(,ff = 30, OOOK, with a 
distance d ~500pc; and assuming Log{g) = 9.5 we found a WD temperature of 32, OOOK with a 
distance of ~300pc. EY Cyg does not show any sign of quasi-molecular satellite lines, and this 
is consistent with our higher temperature model (see the Discussion section), and no need to add 
a second component since this component was driven by the fiux excess which the re-processing 
removed. 



4.2. SS Aurigae 



^SS Aur is a U Gem type DN with an orbital period Porfe=4.3872 hrs (IShafter fc Harkness 

198d ). a WD mass M^a = 1.08 ±0 40 M (7), a secondary mass M2 = 0.39 ± 0.02 Mq, and a system 
inclinat ion i = 38° =1= 16° (Shafter Ill983l ). SS Aur has an HS T FGS parallax measurement of 497 



mas (jHarrison et al. Ill999l ) giving a distance of 201pc. 



Lake &: Sion I (|200ll ) analyzed the lUE archival spectra of SS Aur; their best-fit model photo- 
sphere has Te// = 30, OOOK, Log{g) = 8.0, and solar composition abundances; while the best-fit 
accretion disk model has M,,,^ = 1.0 Mq, i = 41°, and M = 10"^° Mgyr-^ 



Using a newly available Hubble FGS parallax, iLake &: Sion I (j200ll ) showed for the first time 
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that the FUV flux of SS Aur during quiescence was dominated by, if not provided entirely by, a 
hot WD, not an accretion disk. This was something of a milestone because the prevaling view was 
that the underlying degenerate was detected in the FUV only in those dwarf novae which clearly 
revealed the broad Lya profile of a WD, namely U Gem, VW Hyi and WZ Sge. 



More recently ISion et al 



(l2004al ) analyzed the FUSE spectrum of the system. Their best 
fit model is a 27,000K WD contributing 73% of the flux with a hot belt (48,000K) contributing 
the remaining 27% of the flux, with a corresponding distance of 267pc. However, a single WD 
component also gives a good fit with a temperature T^„d=33,000K and a distance of 303pc. They 
also suggested that the absence of the Cm (1175A) absorption line might be a sign of composition 
deficit in carbon. It is important to note here that in these previous studies the spectra of SS Aur 
were not dereddened. 



The STIS spectrum of SS Aur was obtained as part of a snapshot program ([Gansicke et al 



20031 ) and the single STIS spectrum is being analyzed also in ISion et al. I (|2008l ). Here, we ana- 



lyzed the combiiied FU SE+STIS spectrum of SS Aur which we dereddened assuming E(B-V)=0.08 
(jBruch &: Engel II 19941 ). The FUSE spectrum of SS Aur was obtained 28 days after outburst, while 
the STIS spectrum was obtained 34 days after outburst. The fit was carried out with a binning of 
0.2 A in the FUSE spectrum and 0.58 A in the STIS spectrum, consequently the relative weight of 
the FUSE spectrum is about 50% larger than that of the STIS spectrum. 

The best fit model for a ~1.1Mq {Log{g) = 8.71) mass WD and a distance of 200pc gives a 
temperature of 31,000K, solar abundances, a rotational velocity of 400km/s and x^=lA7%. This 
fit is the best when fitting the flux level while keeping the distance and radius of the star fixed. If 
we relax the restriction on the distance, the best fit model with Log{g) = 8.71 gives a temperature 
of 33,000K, a distance of about 240pc and a slightly lower (1-472). The difference in the values 
of the between the 31,000K and the 33,000K models is not significant. However, it is this 
discrepancy between the apparent temperature of the WD and its flux that has led to the adoption 
of second model components in other CVs. However, the WD mass of SS Aur is given with a 
large error (zb0.4MQ). We, therefore, computed more models with Log{g) = 8.31 and 8.93. The 
Log{g) = 8.31 model, for a distance of 200pc gave a temperature of 27,000K, while the lowest 
best fit model for that value of Log{g) gave a temperature of 30,000K and a distance of 254pc. For 
the Log{g) = 8.93 model, the distance of 200pc is achieved with a temperature of 34,000K. This 
model, however, is also the lowest best fit model for that value of Log{g) with the least of 
aU (1.471). 



A model fit to the STIS spectrum alone of SS Aur (which was carried out in (|Sion et al 



20081 )) gave the same results, namely a 34,000K WD assuming Log (g) = 8.8. The only difference 
with our modeling is that the STIS spectrum in ISion et al. I (|2008l ) was NOT normalized to the 
Ff/S'E' spectrum, which explain the slight difference in the value of Log{g). We then ran a model to 
fit the FUSE spectrum alone and found total agreement with the FUSE+STIS spectral fits, namely 
that the best fit model is again a T=34,000K WD, spinning at 400km/sec with solar abundances. 
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for Log{g) = 8.93 and a distance of 200pc. We therefore adopted this model as the best model 
fit for SS Aur, as it is produced within the margin of error of the WD mass and it is the simplest 
(one component) model possible. We do not see any need to add a second FUV component in 
the spectral modeling of SS Aur. The presence of a second (flat and featureless) component to 
the spectrum of SS Aur has not been detected as in the case of VW Hyi, and in the present work 
we restrict our modeling to a single WD component. The Log{g) = 8.93 best fit model with a 
temperature of 34,000K is shown in Figure 3. In Figure 4, we show the same model superposed to 
the FUSE spectrum binned at 0.1 A. 

Usually C II (1066A) and Cm (1175A) are in absorption in a 30,000K WD stellar atmosphere 
while Cm (1175A) and Cm (977A) can show emission originating in a hot second component (e.g. 
thin disk, accretion column, BL). From Figure 4, we see C II (1066A and lOlOA) is in absorption 
while Cm (977A) is in emission. Cm (1175A) on the other hand does not show any sign of 
absorption nor emission (it could be a superposition of absorption and emission). We therefore 
decided to fit the carbon abundances by fitting the C ii lines; achieved by simply assuming solar 
abundances (we note that if we assume a low carbon abundance to fit the Cm 1175 feature, then 
the C II absorption features at lOlOA and 1066A are not fitted) . We do not find evidence that the 
WD is deficient in carbon as is the case for EY Cyg, especially since the Civ (1550A) is in strong 
emission in the STIS snapshot. 

In the present work, we have improved the line identifications and the fit ting in the shorte r 



wavelengths without the need for a second component. The second component (ISion et al. Il2004al ) 
gives too much flux in the short wavelengths. The distance we obtained is in very good agreement 
wi th the distance esti mated from the parallax, which is smaller by 50% than the distance obtained 
in 



Sion et al. I (|2004al ) (a result of our dereddening) . 



4.3. VW Hydri 



VW Hyi is a key system for understanding DNs in general. It is one of the closest (jWarner 



19871 . placed it at 65 pc) and brighte st examples of an SU UMa-type DN. It lies along a line of 
sight with a low interstellar column (jPolidan. Mauche &: Wade lll99Cll . estimated the Hi column 
density to be ~ 6 x 10^^ cm~^), which has permitted study of VW H yi in nearly all wavelength 
ranges, including detection in the usually opaque extreme ultraviolet (IMauche 1119961 ). For these 
reasons, it is one of the best-studied systems. Coherent and quasi-coherent soft X- ray oscillations 



have been detected, and its X-ray luminosity has been found to be unexpectedly low (jBelloni et al 



1991; Mauche et al. 



in the FUV band 



199lh. Its boundary la.y e r (BL) has been show n to emit most of its energy 



Pandel et al. 



(|200l l2005l 1: iGodon fc Sion I (|2005l l. VW Hyi is near the lower 
edge of the CV period gap, with an orbital period of 1.78 hrs and a quiescent optical magnitude 
of 13.8. It is a member of the SU UMa class of DNs, which undergo both normal DN outbursts 
and superoutbursts. The normal outbursts last 1-3 days and occur every 20-30 days, with peak 
visual magnitude of 9.5. The superoutbursts last 10-15 days And occur every 5-6 months, with 
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peak visual magnitude re aching 8.5. The mass of the accreting WD was estimated to be 0.63 Mq 
(jSchoembs Vogt Ill98lh. but more rece ntly a gravitational redshift determination y ielded a larger 



mass M^rf 
1996al lbh. 



O.86M0 (jSion et al. 1119971 ). The inclination of the system is ~ 60" (jHuang et al 



VW Hyi was first observed at UV wavelengths with lUE. Based on observ ations of VW Hyi 



i n qui escence, which showed a broad absorption profile centered on Lyman a, iMateo &: Szkodv 
U9M) argued that UV light from the system was dominated by the WD with Teff = 18, 000ifc2,000 
K (ioi Loq(g) = 8). Much higher S/N spectra were obtained with HST. In particular. ISion et al 



(|l995al ) found the WD temperature to be 22,000K, the first metal abundances, with evidence of 
CNO processing of material in the WD photosphere. They also d erived the first rota tional velocity 



for a WD in a CV using the GHRS G160M spectrum of VW Hyi (jSion et al. Ill995bl ) thus showing 



that weak boundary l ayer emissioii from some CVs could not be explained by their WDs rotating 
near Kepler ian speeds. ISion et al. ] (|i996l . bnoih concluded that the temperature of the WD varied 
by at least 2000K depending on the time since the outburst. They also presented evidence for 
rapidly rotating "accretion belt" in the system. The accretion belt was to be understood physically 
as a region of the W D surface spun up by accreting material with Keplerian velocities. Furthermore, 
Sion et al. I (|l997l ). using the GHRS on HST, measured the gravitational redshift of the WD and 



concluded that the mass of the WD was 0.86^0 32^0 and that the rotation rate of the WD was ~ 
400 km s~^. All of these observations were limited to a wavelength range longward of about 1150 A. 
Howeve r, an 820-184 A sp ectrum of VW Hyi was obtained using the Hopkins Ultraviolet Telescope 



(HUT). iLong et al. I (|l996l ) found that the HUT spectrum was reasonably consistent with a WD 
with a temperature of about 17,000K, but that an improved fit to the data at that time could be 
obtained with a combination of emission from a WD and an accretion disk. 

Higher S/N spectra of VW Hyi in quiescence were obtained more recently with FUSE. With 
a practical wavelength range of 904-1188 A, FUSE is sensitive to a second component, since the 
expected flux from a WD with a temperature of about 20,000 K is very different at 950 and 1100 
A. In addition, the FUSE spectral range and high spectral resolution allows the study of a broad 
range of line transitions. In the present work we only use one FUSE data set (B0700201), which 
presents only a small contribution from the second component. Since the second component is being 
studied elsewhere (Long 2007, private communication) using different FUSE data, sets - E1140109, 
E1140112, E1140113), and it represents o nly a few percent of the flux, we neglect its contribution 



in the present work. In a previous work (jGodon et al. II2004I ) we modeled the FUSE spectrum of 



VW Hyi, including a fast rotating hot accretion belt, and obtained a WD temperature of 23,000K. 
However in that work we did not deredden the spectrum, nor did we include the quasi-molecular 
satellite line in the modeling, and the FUSE spectrum was not combined with the STIS spectrum. 

In the present analysis, we included the quasi-molecular lines in the synthetic spectral modeling, 
to fit the FUSE spectra better in the wavelength range ~ 1060 — IO8OA, and we dereddened the 
spectrum assuming E(B-V)=0.01 (iBruch &: Engel Ill994l ). We are aware that because of the very 
low hydrogen column density one actually expects VW Hyi to have a zero reddening. However, 
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(i) such a small reddening value does not affect the results significantly, and (ii) we want to be 
consistent with the reddening values given in the literature and adopted in this work. The FUSE 
spectrum we used here was obtained 11 days after outburst, while there were three STIS spectra in 
the archive which we combined (weighting them by exposure time) to improve the S/N. Two short 
exposure STIS spectra were obtained 1 and 7 days after outburst; however, they are actually very 
similar to the long exposure STIS spectrum obtained 15 days after outburst. The FUSE and STIS 
spectra do match and fit well together. 

Since the WD mass might be somewhere between ~ 0.6 and 0.9Mq, we computed models for 
two different values of Log{g): 8.0 and 8.5. The model with Log[g) = 8.0 was in better agreement 
with the distance of the system (~ 60pc) and gave a temperature of 22,000K, a projected rotational 
velocity of 400km/s and abundances as specified in Table 5. This model is shown in Figures 5 & 
6. In both figures the FUSE and STIS spectra are binned at 0.1 A, and the relative weight of the 
STIS spectrum is a little above twice that of the FUSE spectrum. Due to the excellent fit of the 
Si, C and N lines, we are confident that these abundances are accurate. While the model does not 
fit the flux (from a second component) in the short wavelengths of FUSE, it accurately follows the 
shape of the continuum and lines at wavelengths A ~ 1050A all the way into the longer wavelengths 
of STIS. 

The inclusion of the quasi-molecular opacity did improve the fit in the FUSE range 1055A- 
1183A, but did not change the fit in the STIS spectrum around ~ 1300 — 1400A range, because of 
the temperature of the model. The WD could be massive with Log{g) ~ 8.5 or larger. 

From this fit, it is evident that a WD alone is not enough to reproduce all the features of the 
spectrum. The bottoms of Lya & Ly/3 do not go to zero and indicate the presence of an additional 
component, which may or may not be linked to the broad emission features. As stated earlier, we 
do not intend in this paper to model the second FUV component of the system as this is being 
done elsewhere. We concentrate here on fitting the main continuum, the absorption lines, the quasi- 
molecular satellite lines and combining the FUSE spectrum with the STIS spectrum in order to 
assess the basic parameters of the WD: temperature, abundances, projected rotational velocity, and 
gravity. Our model neglects possible contribution from a quiescent (and optically thin) accretion 
disk. As an additional test we ran models for the fitting of the STIS spectrum alone and found 
complete agreement with the best fit model for the FUSE+STIS spectrum: a Log{g) = 8.0 WD 
with a temperature T = 22, OOOK and a distance of 60pc, therefore confirming our result for the 
combined FUSE+'SiTlS spectrum. 



4.4. EK Trianguli Australis 



EK TrA is a poorly studied southern DN with an orbital period ju s t abo ve 1.3 hrs. It belongs 
to the SU UMa subclass (optical superhump was confirmed by iHassall I (|l985l )) and has often been 
compared to the well studied DN VW Hyi because of their similarity. However, EK TrA has a much 
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longer outburst period than VW Hyi and can therefore easily be observed during deep quiescence 
for the study of its WD. Despite this, however, the mass of its WD is unknown and its dist ance 
is estimated to be d > 180pc fr om the ri o n-det ection of its secon dary (IGansicke et al. 1119971 ). It 



is important to note that while IWarner I (jl987l ) put it at 200pc, iHassall I (|l985l ) initially derived 
d^ 68 - 86pc. 

EK TrA was observed in August 1980 with lUE and 8 spectra were obtained with the SWP 
and LWR cameras covering mid-outburst and late decline from the outburst. The WD was re- 
vealed and cont ributed about 25% of the lUE SWP flux, with an accretion disk contributing the 
remaining 75% (jGansicke et al. 1119971 ). More recently, EK TrA was observed wit h HST 155 days 



after outburst and a high resolution STIS spectrum was obtai ned and analyzed by (IGansicke et al 



2001 



). They confirmed the temperature of the WD found in IGansicke et al. I () 19971 ) (T=18,800K) 
assuming a canonical mass of O.GMq, derived a rotational velocity of 200km s~^ and slightly sub- 
solar abundances (0.5). They note that assuming a WD mass in the range 0.3-1.4Mq leads to a 
temperature in the range 17, 550-23, 400K. More recently EK TrA was observed with FUSE 45 days 
after outburst and two 33ksec spectra were obtained with the same flux level as the existing STIS 
spectrum. The FUSE spectra have not been modeled previously, and we therefore decided to model 
EK TrA using the combined FUSE+STIS spectrum. 

Since the WD mass of the system is basically unknown, we ran model fits for Log{g) = 7.5, 
8.0, and 8.5 (corresponding to a mass of about 0.40, 0.65 and O.93M0 respectively) and obtained 
a WD temperature of 15,000, 17,000 and 18,000K respectively. The Log{g) = 8.0 best fit model is 
presented in Figure 7. The broad emission lines (in blue/dark grey) have been masked for the fitting 
and both the FUSE and STIS spectra have been binned at 0.175A . In Figure 8 we show the FUSE 
spectrum together with the best fit model shown in Figure 7. In the wavelengths A < 1050Aonly 
the broad O vi doublet and the Ciii(977 A ) emission lines are from the source, all the other sharp 
emission lines are due to air glow. The flux level at these wavelengths is of the same order as 
the noise and we do not attempt to model the spectrum there. In order to fit the line profiles 
in the FUSE range we had to set C=0.1 solar and Si=0.6 solar while keeping S and N to their 
solar values. Here also we found that the inclusion of the quasi-molecular satellite hydrogen line 
opacity improves the fit significantly, especially in the FUSE range where the flux suddenly drops 
shortward of the nitrogen absorption line at 1184 A. In the STIS range this is not as pronounced 
but the drop in flux is also apparent shortward of the Si em ission liries (14 00A). The dis tance we 
obtain ed is at most 137pc {Log{g) = 7.5), somewhere between Hassall ( 19851) estim ate and Warner 



(jl987l ). The temperatu res we obtained here ar e lower than IGansicke et al. I (|l997l ). The differences 



between our result and Gansicke et al 



(jl997l ) are probably due to the fact that, in their analysis, 
Gansicke et al. I (jl997l ) did not include the opacity of the quasi- m olecular satellite hydro gen lines, 
and they did not deredden the spectrum (it has E(B-V)=0.03; iBruch &: Engel I (|l994l )). They 
interpreted the excess of flux in the longer wavelengths of STIS (>1,550A) as possible emission 
from an optically thin accretion disk during quiescence, while we found, after dereddening the 
spectra, no evidence of flux excess. We note here that the complete absence of ISM molecular 
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hydrogen absorption lines in the Ff/5'E' spectrum of EK TrA suggests that the system might indeed 
be located rather nearby (say <150pc). For comparison, in the present work, only the FUSE 
spectrum of VW Hyi at 65pc shows such a complete absence, while SS Aur at 200pc clearly shows 
some molecular hydrogen lines. We therefore speculate that the system is probably located at 
about 125pc, as supported by our results. 



4.5. RU Pegasi 



RU Peg is a U Gem type DN with an orbital period P^rt 
33° , a secondary spectral typ e K2-5V, and a primary (WD) mass M^^i 



8.99 hrs, a system inc lination 

+0.16^ 



1.29" 



1981 



Wade 



1982 



( Stover 



Shaffer Ill983l ). The near-Chandrase khar mass for the WD has been corroborated 



by the Sodium (8190A) doublet radial velocity study of lFriend et al. I (1990). They obtained a mass 
of 1.24 Mq for the WD and also found very good agreement with the range of plausible inclinations, 
found in the studv bv IStover I (Il981h . Recently, a Hubble FGS parallax of 3.55 ± 0.26 mas was 



measured bv I Johnson et al. I (j2003l ) implying a distance of 282pc. 



Sion &: Urban I (j2002l ) analyzed 4 lUE spectra of RU Peg obtained in deep quiescence, and 
found a WD temperature Tg// = 50 — 53, OOOK with a distance of ~ 250 pc, assuming Log{g) = 8.7 
{M = I.IM0). More recently, a FUSE spectrum of RU Peg was obtained 60 days after outburst. 
The FUSE observation lasted for a little more than 3000s (raw exposure time). From all the 5 
systems in this analysis, RU Peg is the most affected by ISM absorption lines. We show the FUSE 
spectrum of RU Peg with line identifications in Figure 9, where we have ma rked only the most 



prominent molecular hydrogen lines. Two modelings have been carried out (jSion et al. Il2004al : 



similar to 



Urban k. Sion I l2006l'l for the F USE spectrum of RU Peg; and both found a WD temperature 



Sion fc Urban I ([2002) assuming M^^ = 1.2Me 



0- 



While we are confident that the WD in RU Peg is actually massive and that its temperature 
is very large (>50,000K), the FUSE spectrum has never been modeled in detail as many of the 
ISM molecular hydrogen lines were misidentified. In addition, since this FUSE analysis was carried 
out, CALFUSE has been updated and we are making use of this implementation. For RU Peg, we 
decided to also model the ISM molecular and atomic hydrogen absorption lines to help improve 
the model fit. 

In Figure 10 we show the best-fit model to the FUSE spectrum of RU Peg (binned at O.lA). 
This model consists of a synthetic WD spectrum model multiplied by the transmission values 
obtained from the ISM model. This WD model has a temperature of 70, OOOK, a rotational velocity 
of 40km/s, Log{g) = 8.8, subsolar carbon and silicon abundances and supra solar nitrogen and 
sulfur abundances as follows. The silicon abundance was set to 0.2 times solar to match the Siiv 
(1066 A) line. The nitrogen abundance was set to solar to match the Niii (980, 1183, & 1184.5 
A) lines. The carbon abundance was set to 0.2 times solar to match the feature at 1108 A and 
the flat portion of the spectrum at 1160-1170 A. The sulfur abundance was set to 10 times solar 
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to match the strong absorption of the S vi and Siv hnes. It is worth noting that if the N/C ratio 
is anomalous in the WD photosphere then this would imply recent thermonuclear processing in 
the WD, bec ause the disk emissio n lines do not indicate a N/C ratio larger than one, but rather 



the opposite (jSion &: Urban II2002I ). However, our abundances determination is not robust because 
the nitrogen and carbon absorption lines are afffected by broad emission lines and ISM molecular 
hydrogen absorption lines. 

The projected rotational velocity was found by matching the profile of the carbon, oxygen, 
silicon and sulfur lines, and in particular the square-shaped bottom of the Cm 1175A absorption 
feature. The ISM model (in Figure 10) has zero metallicity, a temperature of lOOK, a turbulent 
velocity of b=30km s~^, a molecular hydrogen column density of 1 x lO^^cm"^ and an atomic 
hydrogen column density of 1 x lO^^cm"'^. This ISM model fits best the stronger (molecular 
hydrogen) absorption lines but does fit the sharper and weaker (molecular hydrogen) absorption 
lines. 

It is worth noting that all the synthetic WD spectra we generate between T 50, OOOK up to 
T 80, OOOK (with an approximate NLTE line treatment option switched on in SYNSPEC and 
using LTE stellar atmosphere models from TLUSTY) fit the data similarly, the main difference is 
the radius (and therefore mass and gravity) of the WD which is fixed by scaling the model to the 
data assuming a distance of 282pc. Consequently, the 50, OOOK model had Log{g) = 8.55, while the 
70, OOOK model had Log{g) = 8.8. In order to agree with a mass of 1.29M0 we ran models with 
Log{g) = 9.23 with increasing temperatures and found that the model that agrees with the distance 
of 282pc has a temperature of 110, OOOK. At such a high temperature the Siv, Siiv and Oiv lines 
between 1060Aand lOSOAare not deep enough to match the observed features. In addition, since 
we performed only the simplest NLTE approximation, and also obtained \^ = 0.827 (against 0.539 
for the T=70,000K model), we decided not to chose the 110, OOOK model as a best fit model. 

On the other hand, at 50, OOOK, the Oiv line and the nominally strong Niv and Svi lines are 
not pronounced. The temperature had to be increased to 70, OOOK to clearly show the appearance 
of these lines in the synthetic spectrum (including the Oiv line at 1067. SA). Therefore, we chose 
the 70, OOOK model as our best fit model. 

However, the Siiv and Oiv lines (between 1066 and 1068 A) are contaminated by sharp 
ISM absorption lines from molecular hydrogen (L3P2 1066.90A, & L3R3 1067.48A) and the Arl 
1066. 66A line (the other argon line Arl 1048. 20A is also present). However, these sharp ISM lines 
were not modeled accurately in the ISM model presented in Figure 10. Consequently, in order to 
fit the sharper absorption lines of the ISM we generated a second ISM model. This model fits the 
sharp lines better but does not fit the stronger lines as well (it is likely that the ISM has more than 
one component). The second ISM model has zero metallicity, a temperature of 300K, a turbulent 
velocity of b=10km s~^, a molecular hydrogen column density of 1 x lO^^cm"^ and an atomic 
hydrogen column density of 1 x lO^^cm"^ and a blueshift velocity of -35km s^^. In Figure 11 we 
present the 70, OOOK WD synthetic spectrum including the second ISM model in the region IO6OA 
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- 1080A . From this model we see that the O IV hne is distinct from the molecular hydrogen line 
(L3R3) but it is not very strong. As for the Siiv lines, they are contaminated by argon and a 
molecular hydrogen line (L3P2). Since argon (and iron, and all the metals) is not included in the 
ISM model, we cannot confirm whether the Siiv line is present. We also note that Siv is shifted 
to the red compared to the WD model. Indeed, a close look at all the lines from the system reveals 
that some of the (broader) lines are shifted by as much as O.SA to the red; this is especially the 
case for the Svi, Ovi, Siv and Cm lines. This could be due to the motion of the WD in the 
binary, as the red shift is close to the maximum radial velocity shift of the WD in the system and 
the FUSE spectrum had a duration much shorter than the binary period. It is possible that these 
strong absorption lines do not form in the stellar photosphere, but possibly in a corona above it 
or in a hot boundary layer in front of it. This can be determined only if we have line velocities 
in phase-resolved FUV spectra to compare with the velocity derived for photospheric lines in the 
white dwarf's rest frame. 



4.6. The Quasi-Molecular Satellite Lines 



One of the improvements we made in the present work is the modeling of the quasi-molecular 
satellite lines of hydrogen in the FUSE spectra of VW Hyi and EK Tr A. Such model ing has been car- 



ried o ut previously for CV W Ds such as WZ Sge (ISion et al. Ill995al') and AM Her (IGansicke et al 



20061 ). and also for DA WD s (|Koester et al. Ill99d . Il998l : IDuduIs et ahlbooal . bood ) and ZZ Ceti 



WDs (lAllard et al. I l2004bl ) . It is interesting to compare our results with these results as the 



presence/absence of the quasi-molecular hydrogen lines in FUV spectra of WDs is not yet fully 
understood. Quasi-molecular satellite transitions take place during close collisions of the radiating 
hydrogen atom with a perturbing atom or proton. These absorption features are present in the red 
wings of the Lyman series lines (Lya in STIS and lUE spectral ranges at 1400 & I6OOA ; Lyp at 
1058 & 1076A and Ly7 around 995A - in the FUSE spectral range) and provide a noticeable source 



of opa cit y in the hydrogen-ri ch atmosphere of WDs (jKoester et al 



, , , , , , 19851 : lAllard et al. 111994 11998I . 

2004al lbl:l Dupuis et al. 1120061 ). In theory the satellite appearance is very sensitive to the degree 



of ionization and may be used as a temperature diagnostic. On the other hand, since the quasi- 
molecular opacity is proportional to the proton and Hi densities, the quasi-molecular satellite lines 
are expected to be observed at higher effe ctive temperatures in WDs with larger surface gravities 



( Koester et al. 



1996 



Dupuis et al. II2003I ). We summarize in Table 6 the occurrence of the quasi- 



molecular satellite lines in a number of WDs, covering a wide range of temperatures and gravities, 
listed in order of increasing temperature. We are aware that these WDs belong to different types 
(CVs, DA and ZZ Ceti) and that the detection of the quasi-molecular satellite line also depends 
on the S/N of the data. However, this table might help put restrictions on the upper and lower 
limits of the gravity and temperature of WDs, especially when the mass and distance are relatively 
unknown. For example it is clear that the WD of FY Cyg must have Log{g) < 9.3, otherwise one 
would expect to see the quasi-molecular satellite lines in the FUSE range. For WW Cet the oppo- 
site is true and one can therefore expect Log{g) > 7.8. In fact one could argue that for WW Cet 
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Log{g) > 8.3 while for EY Cyg Log{g) < 9.0, as this would be even more consistent with Table 6 
and with the errors on the WD masses of these systems. However, the system that really stands out 



in the list is AM Her which exhibits the 1076A feature but not the one at 1058A . Gansicke et al. 



(j2006l ) suggest that the magnetic field plays a role in the formation of these lines, and cite the fact 
the magn etic DA WD PG 16 58+441 has a weaker than expected 1058A H2 absorption feature. 
However, iDupuis et al. I (|2003l ) mentioned that the treatment of the quasi- molecular satellite opac- 
ities could be improved, especially in the range of ionic densities encountered in the atmosphere 
of ultramassive WDs. The FUSE spectrum of PG 1658+441 actually shows that both the 1058A 
and 1076A absorption features are rather weak, contrary to AM Her where the 1058A feature is 
completely absent. We see the importance of generating a large database of WDs exhibiting the 
quasi- molecular satellite lines of hydrogen, as it could be used to assess the WD's temperature 
and/or gravity more accurately. 



5. Summary 

We have analyzed the FUV spectra of five DNs in quiescence with the aim of securing more 
accurate values of system parameters, such as the WD surface temperatures T^fj, their projected 
rotational velocities Vrotsini (derived mainly from the FUSE data), their chemical abundances, 
their surface gravities and the distances to the systems. Four of these systems (EY CYG, SS Aur, 
VW Hyi, and EK TrA) were part of our broader HST archival program with matching FUSE 
and HST/STIS spectra, while RU Peg has only one FUSE spectrum. We carried out the analysis 
including options for the treatment of the hydrogen molecular satellite lines (for modeling cooler 
WDs) and NLTE atmosphere models (for modeling hotter WDs), dereddening the spectra and 
identifying the ISM molecular absorption lines in the FUSE spectra. The following improvements 
(in methods and results) were achieved. 

For EY Cyg, we disregarded the noisy portion of the FUSE spectrum (A < 950A), which 
turned out to be an artifact of the re-processing and thus improved the fit. From the combined 
HST/STIS + FUSE spectrum we found Tg// = 30,000K and d=530pc assuming Log{g) = 9; and 
Tf.ff = 32, OOOK and d=310pc assuming Log{g) = 9.5. We derived for the first time evidence of 
CNO processing in the photosphere of the WD itself in EY Cyg. Previously, only the emission 
line strengths implied suprasolar N and subsolar (deficient) C. Therefore EY Cyg's WD is only the 
third WD known to have anomalous N/C in its photosphere, the other objects being VW Hyi and 
U Gem. 

For SS Aur, we re-evaluated and improved the WD parameters by modeling for the first time 
for this object the combined FUSE+STIS spectra. We improved the fitting in the FUSE short 



wavelength range without the need to introduce a second component, contrary to (|Sion et al 



2004al )'s result where the sec ond component actually gave too much flux. This might be due to the 



fact that ISion et al. I (l2004al ) did not deredden the FUSE spectrum of SS Aur, while we assumed 



E(B-V)=0.08. In addition this result is in excellent agreement with the distance and we found no 
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deficiency in carbon. We also improved the line identifications over the previously published FUSE 
analysis. 

For VW Hyi, the FUSE and STIS spectra were combined together for the first time. The opac- 
ity of the quasi-molecular satellite lines of hydrogen were included, and the spectrum was dered- 
dened assuming E(B-V)=0.01. Our best model fit has re//=22,000K, Log{g) = 8.0 and d = 60pc. 
Both the projected rotational velocity (400km/s) and the chemical abundances are determined quite 
accurately due to the high S/N of the FUSE spectrum. We find a ratio C/N=0. 25/3=1/12 and the 
silicon abundance is almost twice solar. A second component is detected in the very short wave- 
lengths contributing only a few percent of the total flux. Even though we deliberately decided not 
to model the second componei it, our results are in cornplete agreernent with all the results from pre- 



1996 



Sion et al. 



1996 



1997 



vious analyses of the system (IMateo fc Szkodv Ill984l: ISion et al. 



2001 



Godon et al. 



200- 



1995a 



Gansicke &: Beuermann 



For EK TrA we re-evaluated and improved the WD parameters by modeling for the first time 
the combined FUSE+STIS spectrum. We found a W D temperature o f 17,00 0Kifcl,000K assuming 



L'Og{g) = 8.0 it 0.5 leading to a distance of ~125pc. iGansicke et al. I (|200ll ) carried out a similar 
analysis, but found a somewhat higher temperature and a flux excess in the longer wavelengths, 
assuming a distance d>180pc. The difference in results is probably due to our different modeling; 
namely we combined the FUSE and STIS spectra, included the opacity of the quasi-molecular 
satellite lines for hydrogen, and dereddened the spectrum. 

For RU Peg we improved the line identiflcations over the previously published FUSE analysis 
and we show that the large mass of the WD gives a better agreement with a higher temperature 
than previously estimated; namely the Log{g) = 8.8 model gives a temperature of 70,000K for a 
distance of 282pc. We flnd that the projected rotational velocity is rather small (<70km/s) and 
the carbon and silicon abundances are sub-solar; while sulfur and nitrogen are over- abundant. We 
assessed the molecular hydrogen column density {N{H2) = 1 x lO^^cm"^), the atomic hydrogen 
column density {N{HI) = 1 x lO^'^ - 1 x 10^*^ cm ), and found that the lines are blue shifted with 
a velocity of -35km s~^. 
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Table 1: Parameters of the Systems 



Name 


Type 


Eb~v 


d 




i 


Spectral 


Ml 


Ma 


^ max 


^ min 








(pc) 


(hrs) 


(deg) 


Type 


(Mo) 








EYCyg 


DN UG 


0.00 


450 


11.02377 


16 ± 1 


K7 V 


1.26 ±0.29 


0.59 ±0.13 


11.4 


15.5 


SS Aur 


DN UG 


0.08 


200 


4.3872 


38± 16 


Ml V 


1.08 ± 0.40 


0.39 ±0.02 


10.5 


14.5 


VW Hyi 


DN SU 


0.01 


65 


1.783 


60± 10 


LO 


0.63 - 0.86 


0.11 ±0.02 


8.5 


13.8 


EK TrA 


DN SU 


0.03 


180 


1.5091 


58 ±7 




0.46 ±0.10 


0.09 ± 0.02 


12.0 


17.0 


RU Peg 


DN UG 


0.00 


282 


8.9904 


33 ±5 


K2-5V 


1.29 ±0.20 


0.94 ± 0.04 


9.0 


13.1 



References. EY Cve: Sarna et al. ( 


1995 


; 


Smith et al. ( 199 


7): Costero et al. (199i 


B): Tovmassian et al. 


(2002 


); 


Costero et al. ( 


2004): 


3ion et al. ( 


2004b); SS Aur: 


Shatter 


(1983|); 


Shafter & Harkness 


1986); 


Harrison et al. ( 


19991): VW Hvi: Schoembs & Voet (1981): 


van Amerongen et al. 


1987alb): 


Sion et al. 


(1995b) 


Huang et al. ( 


1996b); 


Long et al. 


( 


19961): Wheatlev et al. 


19961): 


Sion et al. ( 


1997 


); EK TrA: 


VoEt & Semeniuk 


(1980f): Warner ( 


1987): 


Gansicke et al. ( 


1997, 


2001); RU Peg: ^ 


jtover 


(1981): 


Wade 


(1982 


); 


Shatter ( 


198,lh: 


Friend etal. (1990):. 


Johnson et al. 


2003 


). The E(B-V) values are fro 


m 1 Verbunt 


(1987 


); 


La Dous 


(1991) 


Bruch & Engel 


199 


4). 
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Table 2: FUV Observations of the Systems: HST fc FUSE Spectra 



System 


Date 


Telescope 


Exp time 


Dataset 


Filter /Grating 


Operation 


r ^ 




(dd/mm/yy) 


/Instrument 


(ksec) 




/Aperture 


Mode 


(days) 


EYCyg 


16-07-03 


FUSE 


18 


D1450101 


LWRS 


TTAG 


dq*^ 




10-02-03 


STIS 


0.7 


O6LI0V010 


G140L/52X0.2 


ACCUM 


dq^ 


SS Aur 


13-02-02 


FUSE 


14.5 


C1100201 


LWRS 


TTAG 


28 




20-03-03 


STIS 


0.6 


O6LI0F010 


G140L/52X0.2 


ACCUM 


30 


VW Hyi 


18-08-01 


FUSE 


18 


B0700201 


LWRS 


TTAG 


11 




01-06-00 


STIS 


2.5 


O5E205010 


E140M/0.2X0.2 


ACCUM 


1 




06-06-00 


STIS 


2.5 


O5E206010 


E140M/0. 2X0.2 


ACCUM 


7 




10-12-01 


STIS 


9.1 


O5E201010 


E140M/0.2X0.2 


ACCUM 


15 


EK TrA 


24-06-02 


FUSE 


33 


Z9104301 


LWRS 


TTAG 


45 




24-06-02 


FUSE 


33 


Z9104302 


LWRS 


TTAG 


45 




25-07-99 


STIS 


4.3 


O5B612010 


E140M/0.2X0.2 


TTAG 


155 


RU Peg 


05-07-02 


FUSE 


2.9 


CllOOlOl 


LWRS 


TTAG 


60 



"in the last column Tq is the time in quiescence, i.e. since outburst 

''for the very long recurrence time systems (EY Cyg & BZ UMa), the observations were made in deep quiescence (dq) 
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Table 3: FUSE Lines 



Line 

Identification 


Wavelength 

(A) 


Origin 


EY Cyg 


SS Aur 


VW Hyi 


EK TrA 


RU Peg 


Hi 




c,ism 


e,c 


e,c 


e,c 


e,c 


e,c 


Oi 




c,ism 


e,c 


- 


- 


e,c 


- 


Niv 


921.5-924.9 


s 


e,c 


- 


- 


- 


a 


Svi 


933.5 


s 


- 


- 


- 


- 


a 


Svi 


944.5 


s 


- 


- 


- 


- 


a 


Cm 


977.0 


s 


- 


e 


e 


e 


e 


Niii 


979.8 


s 


- 


- 


- 


- 


a 


Nm 


991.6 


s 


c 


c 


- 


- 


c 


Cn 


1010. 


s 


- 


a 


- 


- 


- 


Sin 


1020.7 


ism 


- 


a 


- 


- 


- 


Ovi 


1031.9 


s 


e 


e 


e 


e 


a,e 




1037.6 


s 


e 


c 


e 


e 


a,e 


Cii 


1036.3 


ism 


a 


- 


a 


a 


- 


Siv 


1062.6 


s 


- 


c 


- 


- 


a 


Cii 


1066. 


s 


- 


a 


a 


- 


- 


Si IV 


1066.6 


s 


- 


- 


- 


- 


a 


Oiv 


1067.8 


s 


- 


- 


- 


- 


a 


Siv 


1073.0,1073.5 


s 


a 


a 


- 


- 


a 


Nil 


1084.0-1085.7 


s,c,ism 


c 


c 


a 


c 


c 


Sim 


1108.4-1113.2 


s 






a 


a 




Cm 


1125.6 


s 


a 






a 




Cii 


1127 


s 








a 




Si IV 


1122.5 


s 


a 


a 


a 


a 


a 




1128.3 


s 


a 


a 


a 


a 


a 


Nl 


1134.2-1135.0 


s,c,ism 


c 


c 


c 


c 


c 


Sim 


1140.5-1145.7 


s 


a 


a 


a 






Sim 


1158.1 


s 








a 




Hen 


1168.6 


c 


e,c 




e,c 


e,c 


e,c 


Cm 


1174.9-1176.4 


s 






a 


e,a 


e,a 



The following abbreviations have been used: a - for absorption; e - for emission; c - for contamination (e.g. air glow); 
s - for source; and ism - for interstellar medium. 
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Table 4: STIS Lines 


Line 


Wavelength 


Origin EY Cyg 


bb Aur 


\TWT TJtt^ 

V w iiyi 




Identification 


(A) 












TTT 

Kj III 


1174.9-1176.4 






a 


e 


Dl II 


iiy(j.4 




a,e.' 






^9 

e.' 








a,e? 






^9 

e.' 




iiy4.o 




a,e? 






^9 

e.' 


TTT 

bun 






e 


e 




e 


XJ T 

rl I 


izio. 1 




a,c 


a 


a,e,c 


a,c 


"Mir 
IN V 


Izoo.O 




e 


e 




e 




iz4z.o 




e 


e 




e 


Ol II 


izoy.y, 


iZDU.4 




a 


a 


a 




1zd4.7, 


IZDO.O 




a.' 


a 


a 


T 

Kj I 


iZDD.4 




a.' 




a 




C; TTT 

Dl III 


1 one n 




a 


a 


a 


.^9 


TT 

bi II 


ioU4.4 




a 




a 


r^9 

a.' 


AT T 
IN I 


ioib.U 














1 Q1 Q 7 












T 

Kj I 


lozo.o, 


lozy.O 






a 




TT 
Vj II 


1334.5, 


1335.7 


a 


a,e? 


a 




<^ 1+ 


W1350 








a 


a,e 


C4 T^ ^ 

bi IV 


1393.8 




e 


a,e? 


e,a? 


e 




1402.8 




e 


a,e? 


e,a? 


e 


AT^ TT 

INl II 


1467.3, 


1467.7 










C? T 

b I 


1485.6, 


1487.2 






a 




"NT T 


1492.6, 


1492.8 






a 






1494.7 








a 




Sill 


1526.7 








a 


a 




1533.4 








a 


a 


Civ 


1548.2 






e 


e 


e 




1550.8 






e 


e 


e 


Ci 


1560.3-1561.4 






a 




Hen 


1640.5 




e 






e 


Ci 


1656.3-1658.1 






a? 


a 


Alii 


1670.8 








a? 


a 



The following abbreviations have been used: a - for absorption; e - for emission; c - for contamination (e.g. air glow); 
s - for source; and ism - for interstellar medium. 
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Table 5: Synthetic Stellar Spectra 



Name 


Log{g) 


T 


Vrotsin{i) 


[C] 


[Si] 


[N] 


[S] 


d 


P 


Xi^ 


Figure 




(cm s~^) 


(lO^K) 


(km s~^) 


Solar 


Solar 


Solar 


Solar 


(pc) 








EY Cyg 


9.00 


30.0±2.0 


100±50 


0.03±0.02 


0.6±0.1 


2.0 ±1.0 


6.0±1.0 


530±60 


1.50 


1.140 


1,2 




9.50 


32.0±2.0 


100±50 


0.03±0.02 


0.6±0.1 


2.0 ±1.0 


6.0±1.0 


310±40 


1.50 


1.046 




SS Aur 


8.31 


27.0±1.0 


400±100 


1.0 


1.0 


1.0 


1.0 


200 


0.93 


1.476 






8.31 


30.0±1.0 


400±100 


1.0 


1.0 


1.0 


1.0 


254 


0.93 


1.472 






8.71 


31.0±1.0 


400±100 


1.0 


1.0 


1.0 


1.0 


200 


0.93 


1.474 






8.71 


33.0±1.0 


400±100 


1.0 


1.0 


1.0 


1.0 


240 


0.93 


1.472 






8.93 


34.0±1.0 


400±100 


1.0 


1.0 


1.0 


1.0 


200 


0.93 


1.471 


3,4 


VW Hyi 


8.00 


22.0±1.0 


400±100 


0.25±0.05 


1.8±0.2 


3.0±1.0 


1.0 


60 


1.00 


0.734 


5,6 




8.50 


24.0±1.0 


400±100 


0.25±0.05 


1.8±0.2 


3.0±1.0 


1.0 


51 


1.00 


0.770 




EK TrA 


7.50 


15.5±0.5 


200±100 


0.1±0.05 


0.6±0.1 


1.0 


1.0 


137 


1.00 


0.781 






8.00 


17.0±0.5 


200±100 


0.1±0.05 


0.6±0.1 


1.0 


1.0 


126 


1.00 


0.781 


7,8 




8.50 


18.0±0.5 


200±100 


0.1±0.05 


0.6±0.1 


1.0 


1.0 


104 


1.00 


0.780 




RU Peg 


8.55 


50.0±2.0 


40/-15+30 


0.2±0.1 


0.2±0.1 


>1.0 


10±5 


282 




0.565 






8.72 


60.0±5.0 


40/-15+30 


0.2±0.1 


0.2±0.1 


>1.0 


10±5 


282 




0.545 






8.80 


70.0±5.0 


40/-15+30 


0.2±0.1 


0.2±0.1 


>1.0 


10±5 


282 




0.539 


10,11 




8.92 


80.0±5.0 


40/-15+30 


0.2±0.1 


0.2±0.1 


>1.0 


10±5 


282 




0.553 






9.23 


110.0±5.0 


40/-15+30 


0.2±0.1 


0.2±0.1 


>1.0 


10±5 


282 




0.827 





"The factor by which the STIS spectrum had to be scaled so that its flux level (continuum) matches the flux level 
(continuum) of the corresponding FUSE spectrum at the overlapping wavelengths w 1150 — 1185A 
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Table 6: Presence of Quasi-Molecular Satellite Lines {H2) in White Dwarfs 



Name 


type 


Teff 


Log{g) 


FUSE range 


STIS range 






(lOOOK) 




1058,10761 


1,400A 


G226-29 


ZZ Ceti 


12.0 


7.93 


yes, yes 


_ 


G185-32 


ZZ Ceti 


12.0 


7.9 


yes, yes 


- 


WZ Sge 


DN UG 


14.0 


8.5 




yes 


WZ Sge 


DN UG 


15.5 


8.5 




weak 


40 Eri B 


DA WD 


16.5 


7.77 


yes, yes 




EK TrA 


DN SU 


17.0 


8.0 


yes, yes 


weak 


BZ UMa 


DN SU 


17.5 


8.5 




yes 


AM Her 


NL AM 


19.8 


8.2 


no, yes 


very weak 


Wolf 1346 


DA WD 


19.9 


7.84 


yes, yes 




He 3 


DA WD 


21.4 


8.04 


yes, yes 




GD 140 


DA WD 


21.6 


8.41 


yes, yes 




VW Hyi 


DN SU 


22.0 


8.0 


yes, yes 


no 


WZ Sge 


DN UG 


23.2 


8.5 


yes, yes 




L825-14 


DA WD 


25.0 


7.76 


no, no 




WW Get 


DN UG 


26.0 


8.3 


yes, yes 


no 


PG 1658+441 


DA WD 


29.6 


9.31 


weak, yes 




EYCyg 


DN UG 


30.0 


9.0 


no, no 


no 



References. G226-29 & G185-32: Ullard et aTl (|2004bl ): BZ UMa: based on an inspe ction of the STIS spectrum com- 



pared to a synthetic spectrum; EK TrA: this work; AM Her: iGansicke et al 
the quasi-m olecular feature was detected but not modeled in 
( 200: 



Godon et al. 



20061): all the other DA WDs are from the ORFEUS observations of 



(2006 



(2006 



VW Hyi: this work; WW Get 
PG 1685+441: IPupuis et aF 



Koester et al. ] (|l998l) . WZ Sge: all the 



modelings carried out in ISion et al. I (|l995ah : iGodon et all (|2006al ): iLong et al. I \200i ) assumed a 0.9 solar mass 
The 2004 July STIS spectrum of WZ Sge (with Te// = 15,500K) was analyzed again in the present work to detect the 
quasi-molecular hydrogen feature around 1400A : it was found to be rather weak. 
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FIGURE CAPTIONS 
Fig. 1.— 

The combined FUSE+STIS spectrum of EY Cyg (light grey/red Une) is shown together with a 
30,000K single WD model (thick black line) assuming Log{g) = 9.0. The regions masked during 
fitting are shown in blue (dark grey). With a flux of only a few 10~^^ergs s~^cm~^A~^, EY Cyg 
is a relatively weak source contaminated with many air glow emission lines in the FUSE range, 
especially in the shorter wavelength. The increase of flux at A < 934A is masked in the modeling 
(see text). The WD model has a projected rotational velocity of 100km s~^, a very low carbon 
abundance (note the absence of the Cm (977A) and CiV (1550A) emission lines ) and xt = 1-108. 
The distance obtained from fitting the fluxes is ^ 500pc. 

Fig. 2.— The 30,000K WD model (from Fig.l) is shown here with the FUSE spectra of EY Cyg 
(binned here at 0.1 A ) and line identifications. Air glow lines are annotated with a plus sign inside a 
circle; the O i, N i and N II emission lines are also due to air glow; the interstellar H2 molecular lines 

have been labeled vertically; the O VI doublet emission is from the source; - all these emission lines 
together with the ISM molecular hydrogen absorption lines have been masked in the fitting. The 
abundance of C, S, N and Si have been set to 0.03, 6.0, 2.0 and 0.6 solar (respectively) by fitting 
the following absorption lines: SiV (1173A), Cli (lOlOA, 1066A), Cm (1175A), Sim (1140-1145A) 
and Nil (1184-1186 A). 

Fig. 3. — The combined FUSE+STIS spectrum of SS Aur (in red/light grey) together with the 
best fit WD model (in black). The regions that have been masked for the fitting are shown in blue 
(dark grey). The synthetic stellar spectrum has a temperature of 34,000K (assuming Log{g) = 8.93 
and a distance of 200pc), a projected rotational velocity of 400km and solar abundances. The 
spectrum has been dereddened assuming E(B-V)=0.08. 

Fig. 4. — The FUSE spectrum of SS Aur (binned here at 0.1 A ) is shown with the best fit model 

shown in Figure 3. Note that the C iii (977 A) emission indicates that there could possibly be some 
emission at 1175A, which would make it difficult to detect absorption there, potentially explaining 
the observed flat spectrum in that region. For that reason the carbon abundance was fltted using 
the C II lines seen in absorption both at 1010 A and 1066 A. The sharp emission lines due to air 
glow and the ISM molecular hydrogen absorption lines have been masked for the fitting. 

Fig. 5. — The combined FUSE+STIS spectrum of VW Hyi is shown in red (light grey) together with 
the best fit WD model (in black). The emission lines that have been masked are shown in blue (dark 
grey). The synthetic stellar spectrum has a temperature T = 22, OOOif (assuming Log{g) = 8.0), a 
projected rotational velocity of Vrotsin{i) = 400km s~^, and non-solar composition: C=0.25x solar, 
Si=1.8xsolar, N=3.0xsolar. The spectrum has been dereddened assuming E(B-V)=0.01 and the 
distance obtained from the fitting is 60pc. 

Fig. 6. — The FUSE portion of Figure 6 is shown with line identifications. The dotted-dashed 
line shows the same model without the inclusion of the quasi-molecular satellite lines at 1058A and 
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1076A (shown with arrows). The synthetic spectrum does not attempt to model the O vi doublet 
and C III (977A) emission lines and the low-level flux flat continuum in the lower wavelengths. 

Fig. 7. — The combined FUSE+STIS spectrum of EK TrA is shown in red (light grey) together 
with the best fit WD model (in black). The strong emission lines that have been masked for 

the fitting are shown in blue (dark grey). The synthetic stellar spectrum has a temperature T = 
17, OOOK (assuming Log{g) = 8.0), a projected rotational velocity of 200km s~^, and the following 
abundances (in solar units): C=0.1 and Si=0.6. The spectrum was dereddened assuming E(B- 
V)=0.03 and the distance obtained from the fitting is d=126pc. 

Fig. 8. — The FUSE spectrum of EK Tra (in red/light grey) is shown with the best fit model 
presented in Figure 6 and line identifications. The dotted-dashed line shows the same model 
without the inclusion of the quasi-molecular satellite lines modeling, at 1058A and 1076A . With 
a flux of «i 1 X 10""-'^^ergs s^^cm^^A"-*^ and less, EK TrA is actually a weak FUSE source. All 
the sharp emission lines arc due to contamination (helio- and geo-coronal emission) and have been 
masked for the fitting. The broad Cm (977A & 1175A ) and Ovi (~ 1130- 1140A ) emission lines 
arc associated with the source. The synthetic spectrum does not attempt to model the O VI doublet 
and C III emission lines and the low-level flux flat continuum in the lower wavelengths (which is of 
the level of the noise). The Cll (1036A) absorption line is interstellar. 

Fig. 9. — The FUSE spectrum of RU Peg in quiescence. The sharp ISM molecular hydrogen 
absorption lines have been annotated as in Figure 2. The very broad emission lines of the Ovi 
doublet and Cm (977 & 1175 A) are clearly apparent. In addition there are absorption lines which 
originate only at very high temperatures (T > 50, OOOK) such as Niv (around 923A), S VI (933.5 & 
944. 5A), and Ovi (1031.7 & 1037.3A). The additional sharp absorption lines (from low ionization 
levels such as Nl) are either from the ISM or due to air contamination (or both). Though the flux 
is « 1 X 10~^^ergs s~^cm~^ A~^, the spectrum is clearly under-exposed (Iksec) and has a low S/N. 

Fig. 10. — The best fit WD model to the FUSE spectrum of RU Peg is shown in this Figure. 
The synthetic spectrum is in solid black, the observed spectrum is in red (light grey), the regions 
of the observed spectrum masked for the fitting are shown in blue (dark grey). The synthetic 
stellar spectrum consists of 70, OOOK WD, with a projected rotational velocity of 40km s^^, carbon 
abundance of 0.2 times solar, silicon abundance of 0.2 times solar, and sulfur abundance of 10 times 
solar. The absorption lines that have been used to produce the abundance fits are annotated. 

Fig. 11. — A detail of the fitting of the absorption lines around 1070A (note the Y-axis scale!). The 
molecular hydrogen lines have been marked below the spectrum for clarity. The absorption lines 
shown are SiV (1062. 7A), Si IV (1066.6A), OlV (1067.8A), and the SiV doublet (around 1073.5A). 
The ISM model does not include metals and therefore the An and Fell lines are not modeled, 
though they are present in the observed spectrum. Note the redshift of the SiV doublet. 
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